Converging and evolving immuno-genomic routes toward immune escape in breast cancer.

The interactions between tumor and immune cells along the course of breast cancer progression remain largely unknown. Here, we extensively characterize multiple sequential and parallel multiregion tumor and blood specimens of an index patient and a cohort of metastatic triple-negative breast cancers. We demonstrate that a continuous increase in tumor genomic heterogeneity and distinct molecular clocks correlated with resistance to treatment, eventually allowing tumors to escape from immune control. TCR repertoire loses diversity over time, leading to convergent evolution as breast cancer progresses. Although mixed populations of effector memory and cytotoxic single T cells coexist in the peripheral blood, defects in the antigen presentation machinery coupled with subdued T cell recruitment into metastases are observed, indicating a potent immune avoidance microenvironment not compatible with an effective antitumor response in lethal metastatic disease. Our results demonstrate that the immune responses against cancer are not static, but rather follow dynamic processes that match cancer genomic progression, illustrating the complex nature of tumor and immune cell interactions.

Genomic and epigenomic basis of breast invasive lobular carcinomas lacking CDH1 genetic alterations.

CDH1 (E-cadherin) bi-allelic inactivation is the hallmark alteration of breast invasive lobular carcinoma (ILC), resulting in its discohesive phenotype. A subset of ILCs, however, lack CDH1 genetic/epigenetic inactivation, and their genetic underpinning is unknown. Through clinical targeted sequencing data reanalysis of 364 primary ILCs, we identified 25 ILCs lacking CDH1 bi-allelic genetic alterations. CDH1 promoter methylation was frequent (63%) in these cases. Targeted sequencing reanalysis revealed 3 ILCs harboring AXIN2 deleterious fusions (n = 2) or loss-of-function mutation (n = 1). Whole-genome sequencing of 3 cases lacking bi-allelic CDH1 genetic/epigenetic inactivation confirmed the AXIN2 mutation and no other cell-cell adhesion genetic alterations but revealed a new CTNND1 (p120) deleterious fusion. AXIN2 knock-out in MCF7 cells resulted in lobular-like features, including increased cellular migration and resistance to anoikis. Taken together, ILCs lacking CDH1 genetic/epigenetic alterations are driven by inactivating alterations in other cell adhesion genes (CTNND1 or AXIN2), endorsing a convergent phenotype in ILC.

Nonlobular Invasive Breast Carcinomas with Biallelic Pathogenic CDH1 Somatic Alterations: A Histologic, Immunophenotypic, and Genomic Characterization.

CDH1 encodes for E-cadherin, and its loss of function is the hallmark of invasive lobular carcinoma (ILC). Albeit vanishingly rare, biallelic CDH1 alterations may be found in nonlobular breast carcinomas (NL-BCs). We sought to determine the clinicopathologic characteristics and repertoire of genetic alterations of NL-BCs harboring CDH1 biallelic genetic alterations. Analysis of 5842 breast cancers (BCs) subjected to clinical tumor-normal sequencing with an FDA-cleared multigene panel was conducted to identify BCs with biallelic CDH1 pathogenic/likely pathogenic somatic mutations lacking lobular features. The genomic profiles of NL-BCs with CDH1 biallelic genetic alterations were compared with those of ILCs and invasive ductal carcinomas (IDCs), matched by clinicopathologic characteristics. Of the 896 CDH1-altered BCs, 889 samples were excluded based on the diagnosis of invasive mixed ductal/lobular carcinoma or ILC or the detection of monoallelic CDH1 alterations. Only 7 of the 5842 (0.11%) BCs harbored biallelic CDH1 alterations and lacked lobular features. Of these, 4/7 (57%) cases were ER-positive/HER2-negative, 1/7 (14%) was ER-positive/HER2-positive, and 2/7 (29%) were ER-negative/HER2-negative. In total, 5/7 (71%) were of Nottingham grade 2, and 2/7 (29%) were of grade 3. The NL-BCs with CDH1 biallelic genetic alterations included a mucinous carcinoma (n = 1), IDCs with focal nested growth (n = 2), IDC with solid papillary (n = 1) or apocrine (n = 2) features, and an IDC of no special type (NST; n = 1). E-cadherin expression, as detected by immunohistochemistry, was absent (3/5) or aberrant (discontinuous membranous/cytoplasmic/granular; 2/5). However, NL-BCs with CDH1 biallelic genetic alterations displayed recurrent genetic alterations, including TP53, PIK3CA (57%, 4/7; each), FGFR1, and NCOR1 (28%, 2/7, each) alterations. Compared with CDH1 wild-type IDC-NSTs, NL-BCs less frequently harbored GATA3 mutations (0% vs 47%, P = .03), but no significant differences were detected when compared with matched ILCs. Therefore, NL-BCs with CDH1 biallelic genetic alterations are vanishingly rare, predominantly comprise IDCs with special histologic features, and have genomic features akin to luminal B ER-positive BCs.

KIT genetic alterations in breast cancer.

AIMS: Activating somatic mutations or gene amplification of KIT result in constitutive activation of its receptor tyrosine kinase, which is targetable in various solid tumours. Here, we sought to investigate the presence of KIT genetic alterations in breast cancer (BC) and characterise the histological and genomic features of these tumours. METHODS: A retrospective analysis of 5,575 BCs previously subjected to targeted sequencing using the FDA-authorised Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Targets (MSK-IMPACT) assay was performed to identify BCs with KIT alterations. A histological assessment of KIT-altered BCs was conducted, and their repertoire of genetic alterations was compared with that of BCs lacking KIT genetic alterations, matched for age, histological type, oestrogen receptor/HER2 status and sample type. RESULTS: We identified 18 BCs (0.32%), including 9 primary and 9 metastatic BCs, with oncogenic/likely oncogenic genetic alterations affecting KIT, including activating somatic mutations (n=4) or gene amplification (n=14). All KIT-altered BCs were of high histological grade, although no distinctive histological features were observed. When compared with BCs lacking KIT genetic alterations, no distinctive genetic features were identified. In two metastatic KIT-altered BCs in which the matched primary BC had also been analysed by MSK-IMPACT, the KIT mutations were found to be restricted to the metastatic samples, suggesting that they were late events in the evolution of these cancers. CONCLUSIONS: KIT genetic alterations are vanishingly rare in BC. KIT-altered BCs are of high grade but lack distinctive histological features. Genetic alterations in KIT might be late events in the evolution and/or progression of BC.

Cancer-Causative Mutations Occurring in Early Embryogenesis.

Mosaic mutations in normal tissues can occur early in embryogenesis and be associated with hereditary cancer syndromes when affecting cancer susceptibility genes (CSG). Their contribution to apparently sporadic cancers is currently unknown. Analysis of paired tumor/blood sequencing data of 35,310 patients with cancer revealed 36 pathogenic mosaic variants affecting CSGs, most of which were not detected by prior clinical genetic testing. These CSG mosaic variants were consistently detected at varying variant allelic fractions in microdissected normal tissues (n = 48) from distinct embryonic lineages in all individuals tested, indicating their early embryonic origin, likely prior to gastrulation, and likely asymmetrical propagation. Tumor-specific biallelic inactivation of the CSG affected by a mosaic variant was observed in 91.7% (33/36) of cases, and tumors displayed the hallmark pathologic and/or genomic features of inactivation of the respective CSGs, establishing a causal link between CSG mosaic variants arising in early embryogenesis and the development of apparently sporadic cancers. SIGNIFICANCE: Here, we demonstrate that mosaic variants in CSGs arising in early embryogenesis contribute to the oncogenesis of seemingly sporadic cancers. These variants can be systematically detected through the analysis of tumor/normal sequencing data, and their detection may affect therapeutic decisions as well as prophylactic measures for patients and their offspring. See related commentary by Liggett and Sankaran, p. 889. This article is highlighted in the In This Issue feature, p. 873.

Morphologic and Genomic Characteristics of Breast Cancers Occurring in Individuals with Lynch Syndrome.

PURPOSE: Lynch syndrome is defined by germline pathogenic mutations involving DNA mismatch repair (MMR) genes and linked with the development of MMR-deficient colon and endometrial cancers. Whether breast cancers developing in the context of Lynch syndrome are causally related to MMR deficiency (MMRd), remains controversial. Thus, we explored the morphologic and genomic characteristics of breast cancers occurring in Lynch syndrome individuals. EXPERIMENTAL DESIGN: A retrospective analysis of 20,110 patients with cancer who underwent multigene panel genetic testing was performed to identify individuals with a likely pathogenic/pathogenic germline variant in MLH1, MSH2, MSH6, or PMS2 who developed breast cancers. The histologic characteristics and IHC assessment of breast cancers for MMR proteins and programmed death-ligand 1 (PD-L1) expression were assessed on cases with available materials. DNA samples from paired tumors and blood were sequenced with Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (≥468 key cancer genes). Microsatellite instability (MSI) status was assessed utilizing MSISensor. Mutational signatures were defined using SigMA. RESULTS: A total of 272 individuals with Lynch syndrome were identified, 13 (5%) of whom had primary breast cancers. The majority of breast cancers (92%) were hormone receptor-positive tumors. Five (42%) of 12 breast cancers displayed loss of MMR proteins by IHC. Four (36%) of 11 breast cancers subjected to tumor-normal sequencing showed dominant MSI mutational signatures, high tumor mutational burden, and indeterminate (27%) or high MSISensor scores (9%). One patient with metastatic MMRd breast cancer received anti-PD1 therapy and achieved a robust and durable response. CONCLUSIONS: A subset of breast cancers developing in individuals with Lynch syndrome are etiologically linked to MMRd and may benefit from anti-PD1/PD-L1 immunotherapy.

TERT promoter hotspot mutations and gene amplification in metaplastic breast cancer.

Metaplastic breast cancers (MBCs) are characterized by complex genomes, which seem to vary according to their histologic subtype. TERT promoter hotspot mutations and gene amplification are rare in common forms of breast cancer, but present in a subset of phyllodes tumors. Here, we sought to determine the frequency of genetic alterations affecting TERT in a cohort of 60 MBCs with distinct predominant metaplastic components (squamous, 23%; spindle, 27%; osseous, 8%; chondroid, 42%), and to compare the repertoire of genetic alterations of MBCs according to the presence of TERT promoter hotspot mutations or gene amplification. Forty-four MBCs were subjected to: whole-exome sequencing (WES; n = 27) or targeted sequencing of 341-468 cancer-related genes (n = 17); 16 MBCs were subjected to Sanger sequencing of the TERT promoter, TP53 and selected exons of PIK3CA, HRAS, and BRAF. TERT promoter hotspot mutations (n = 9) and TERT gene amplification (n = 1) were found in 10 of the 60 MBCs analyzed, respectively. These TERT alterations were less frequently found in MBCs with predominant chondroid differentiation than in other MBC subtypes (p = 0.01, Fisher's exact test) and were mutually exclusive with TP53 mutations (p < 0.001, CoMEt). In addition, a comparative analysis of the MBCs subjected to WES or targeted cancer gene sequencing (n = 44) revealed that MBCs harboring TERT promoter hotspot mutations or gene amplification (n = 6) more frequently harbored PIK3CA than TERT wild-type MBCs (n = 38; p = 0.001; Fisher's exact test). In conclusion, TERT somatic genetic alterations are found in a subset of TP53 wild-type MBCs with squamous/spindle differentiation, highlighting the genetic diversity of these cancers.

Genetic characterisation of adult primary pleomorphic uterine rhabdomyosarcoma and comparison with uterine carcinosarcoma.

AIMS: To characterise the genetic alterations in adult primary uterine rhabdomyosarcomas (uRMSs) and to investigate whether these tumours are genetically distinct from uterine carcinosarcomas (UCSs). METHODS AND RESULTS: Three tumours originally diagnosed as primary adult pleomorphic uRMS were subjected to massively parallel sequencing targeting 468 cancer-related genes and RNA-sequencing. Mutational profiles were compared with those of UCSs (n = 57) obtained from The Cancer Genome Atlas. Sequencing data analyses were performed using validated bioinformatic approaches. Pathogenic TP53 mutations and high levels of genomic instability were detected in the three cases. uRMS1 harboured a likely pathogenic YTHDF2-FOXR1 fusion. uRMS2 harboured a PPP2R1A hotspot mutation and amplification of multiple genes, including WHSC1L1, FGFR1, MDM2, and CCNE1, whereas uRMS3 harboured an FBXW7 hotspot mutation and an ANKRD11 homozygous deletion. Hierarchical clustering of somatic mutations and copy number alterations revealed that these tumours initially diagnosed as pleomorphic uRMSs and UCSs were similar. Subsequent comprehensive pathological re-review of the three uRMSs revealed previously unidentified minute pan-cytokeratin-positive atypical glands in one case (uRMS3), favouring its reclassification as UCS with extensive rhabdomyosarcomatous overgrowth. CONCLUSIONS: Adult pleomorphic uRMSs harbour TP53 mutations and high levels of copy number alterations. Our findings underscore the challenge in discriminating between uRMS and UCS with rhabdomyosarcomatous differentiation.

Multilevel biological characterization of exomic variants at the protein level significantly improves the identification of their deleterious effects.

MOTIVATION: There are now many predictors capable of identifying the likely phenotypic effects of single nucleotide variants (SNVs) or short in-frame Insertions or Deletions (INDELs) on the increasing amount of genome sequence data. Most of these predictors focus on SNVs and use a combination of features related to sequence conservation, biophysical, and/or structural properties to link the observed variant to either neutral or disease phenotype. Despite notable successes, the mapping between genetic variants and their phenotypic effects is riddled with levels of complexity that are not yet fully understood and that are often not taken into account in the predictions, despite their promise of significantly improving the prediction of deleterious mutants. RESULTS: We present DEOGEN, a novel variant effect predictor that can handle both missense SNVs and in-frame INDELs. By integrating information from different biological scales and mimicking the complex mixture of effects that lead from the variant to the phenotype, we obtain significant improvements in the variant-effect prediction results. Next to the typical variant-oriented features based on the evolutionary conservation of the mutated positions, we added a collection of protein-oriented features that are based on functional aspects of the gene affected. We cross-validated DEOGEN on 36 825 polymorphisms, 20 821 deleterious SNVs, and 1038 INDELs from SwissProt. The multilevel contextualization of each (variant, protein) pair in DEOGEN provides a 10% improvement of MCC with respect to current state-of-the-art tools. AVAILABILITY AND IMPLEMENTATION: The software and the data presented here is publicly available at http://ibsquare.be/deogen CONTACT: : wvranken@vub.ac.be SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DIDA: A curated and annotated digenic diseases database.

DIDA (DIgenic diseases DAtabase) is a novel database that provides for the first time detailed information on genes and associated genetic variants involved in digenic diseases, the simplest form of oligogenic inheritance. The database is accessible via http://dida.ibsquare.be and currently includes 213 digenic combinations involved in 44 different digenic diseases. These combinations are composed of 364 distinct variants, which are distributed over 136 distinct genes. The web interface provides browsing and search functionalities, as well as documentation and help pages, general database statistics and references to the original publications from which the data have been collected. The possibility to submit novel digenic data to DIDA is also provided. Creating this new repository was essential as current databases do not allow one to retrieve detailed records regarding digenic combinations. Genes, variants, diseases and digenic combinations in DIDA are annotated with manually curated information and information mined from other online resources. Next to providing a unique resource for the development of new analysis methods, DIDA gives clinical and molecular geneticists a tool to find the most comprehensive information on the digenic nature of their diseases of interest.